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Abstract 

It is possible to numerically solve QCD in the planar limit using standard numerical techniques 
on existing computer clusters. The basic ideas behind the computational strategy and recent 
numerical results in the fermionic sector of large QCD are reviewed. 
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I. INTRODUCTION 



The large limit of QCD has several attractive features: 

• Only planar diagrams contribute at fixed g'^N 

• Fermions in the fundamental representation are autornatically quenched as long as the 
number of flavors is kept flnite when taking N ^ oo llj] . 

• It is expected to preserve the essential perturbative and non-perturbative features of 
the iV = 3 theory. The exact solution of the two dimensional theory reveals a tower 
of mesons . 



Phenomenological predictions from large N QCD tend to be close to experimental re- 
sults indicating that leading corrections in the expansion are at least occasionally 
small 



a. 



• Recent work indicates a duality between large N gauge theory and some string theory 
at zero string coupling 

Lattice QCD provides a non-perturbative method to numerically solve large QCD. 
Conventional methods involve computing plr|^sical quantities as a function of and extrap- 
olating these results to the N ^ oo limit 5]. So far, this approach has only been used to 
compute observables in the pure gauge theory, but is has provided evidence that large N 
can be quantitatively useful even at = 3. Conventional computations proceed as follows: 
One flrst takes the large volume limit followed by the continuum limit at a flxed A^ and then 
extrapolates to A^ — oo. This approach will be difficult for fermionic quantities. This is due 
to the fact that the quenched theory at flnite A^ has pathological singularities [gI . Therefore, 
it is necessary to flrst compute fermionic observables in a full dynamical simulation at a 
flxed A^ and then take the N ^ oo limit. 

A short cut to obtain results in the A^ ^ oo limit using lattice QCD should use two 
aspects of large A^ QCD. One is that the physical volume corrections are expected to be 
small in the N —>■ oo limit since one has inflnite number of degrees of freedom at a single 
lattice site. Secondly, fermions are naturally quenched in the N oo limit. To make use 
of the second feature, it is best to take the N ^ oo limit of a quenched theory at a flxed 
lattice coupling, b = l/{g'^N), before taking the large volume limit and the continuum limit. 



Since we are in a fixed lattice volume and lattice coupling, all fermionic observables are 
well-defined and therefore the N —>■ oo limit will be well-defined. Since we are taking the 
large volume limit in the N ^ oo theory, this is also well-defined and does not suffer from 
quenched pathologies. Lastly, we are taking the continuum limit of a set of well-defined 
theories on the lattice. 

The procedure described in the previous paragraph will work independent of the size of 
the finite volume corrections in the — > oo limit. But, as mentioned before, this correction 
is expected to be small. In fact, tremendous excitement was generated by a paper of Eguchi 
and Kawai in the early 1980s since it showed that the lattice theory in the N oc limit 
reduced to a single site model if the U'^{1) symmetries specific to this lattice model remain 
unbroken. This theory has only d SU(N) matrices with N —>■ oo and the [/"^(l) symmetries 
are the N —>■ oo limit of the Zff symmetries associated with the d Polyakov loops in the 
single site model. Unfortunately, the U'^{1) symmetries are broken in the weak coupling 

a single site model is very attractive, solutions 
to fix the above problem were immediately found. One is referred to as the quenched 
Eguchi-Kawai model and another is referred to as the twisted Eguchi-Kawai model Qj. 
The twisted Eguchi-Kawai model is an elegant solution to the problem encountered by the 
original proposal of Eguchi and Kawai and has also led to the discovery of non- commutative 
gauge field theories jlO]. But, it is possible to include fermions only under some very specific 
conditions in the twisted Eguchi-Kawai model. 

Quenched Eguchi-Kawai solves the problem of the original Eguchi-Kawai in an intuitive 
fashion and has other implications. The idea is to restrict the integral over all SU(N) 
matrices with a fixed set of N eigenvalues. Observables are first computed under this 
restriction and these observables are then averaged over all choices for the eigenvalues. 
Freezing the eigenvalues of the SU(N) matrix in the computation of physical observables can 
be interpreted as quenched momentum prescription for gluons shown to be valid in the realm 
of perturbation theo rv In fact, the interpretation of quenched momentum prescription 
extends to fermions [12]. Fermion propagators are computed in the pure gauge theory owing 
to the fact that fermion loops are quenched in the N —>■ oo limit. But, fermion propagators 
with momenta that differ by in a fixed gauge field background are computed by using 
two different gauge field backgrounds that differ by a U{1) phase of e*^'' in the d directions. 
This is referred to as the quenched momentum prescription for fermions. Since we integrate 
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over all SU(N) fields, momentum for a single propagator has no meaning. Only physical 
observables like meson propagators carry momentum where the two quarks that make up 
the meson see gauge fields that differ by the U(l) phase. Furthermore, one sees that mesons 
can have momentum as small as 27r/A^ on a single site lattice and therefore one can study 
mesons for a continuum of momenta on a single site lattice in the N oo limit. 

The proposal to numerically solve QCD in the large N limit should therefore proceed 
as follows. Consider the single site SU(N) lattice gauge theory with the Wilson plaquette 
action. This model has d SU(N) matrices denoted by U^, /i = 1, - ■ ■ ,d. Generate gauge 
fields U according to the single site Wilson plaquette action with the added restriction that 
the d SU(N) integrals go over all SU(N) matrices with a fixed set of eigenvalues that are 
uniformly distributed on the unit circle. Compute quark propagators in the fixed gauge field 
background e^^^^U^ where is restricted to take values equal to such that modified 
gauge field seen by the quark belongs to SU(N). This quark propagator, G{p,U) carries 
a momentum equal to p^. A meson propagator is computed using TrG(p, U)G{q, U) and 
the result is a meson carrying a momentum P = p — q. Averaging over all gauge field 
configurations results in the meson propagator in the large A^ limit of QCD. To recover 
"translational invariance" , it is necessary to sum over all p and q for a fixed P = p — q. This 
amounts to a computation of A^'^ propagators in a fixed gauge field background and this is 
a time consuming numerical computation for d = 4. But the idea of momentum quenching 
can be tested in d = 2 and the results fU] agree with the exact solution of 't Hooft P]. 

Eigenvalues of f/^ are uniformly distributed on the unit circle in the quenched Eguchi- 
Kawai model. Then, e*^f't/^ with = is just a cyclic permutation of all the eigenvalues 
of Ufj^ and should be included in the average over all gauge fields. But the numerical algorithm 
that generates the gauge fields does not succeed in doing this and it is therefore necessary 
to sum over all p and q for a fixed P = p — q when computing meson propagators. Since it 
is necessary to take the large A^ limit first, the proposal in the previous paragraph is quite 
a bit more numerical involved than one naively envisions. Therefore, it calls for a different 
realization of the original Eguchi-Kawai idea. The reduction argument of Eguchi and Kawai 
also applies to the lattice model on an periodic lattice. Reduction will go through and 
there will not be any finite volume effects if the d U{1) symmetries associated with the d 
Polyakov loops are not broken. Indeed, numerical analysis shows that the above symmetries 
are not broken as long as L > Lc{b) for d = 3 and d = 4 Furthermore, Lc{b) 



4 



scales properly as one takes the continuum limit indicating the existence of a physical sea 
Ic such that large QCD is properly reproduced, without finite volume effects if I > L |l5| 
This is referred to as continuum reduction and it results in a modification of the proposal 
in the previous paragraph. One works on an L'^ lattice with L just above Lc{b) instead of 
working on a 1*^ lattice. The gauge fields are generated using the standard single plaquette 
Wilson action with no restriction on the gauge fields themselves. All other aspects of the 
proposal remains the same, with the fermion momenta now going in steps of This 
modified proposal has been shown to work in practice resulting in the computation in the 
chiral condensate in the large N limit of QCD Currently, work is close to completion 



in the computation of the pion decay constant, F^, in the large limit of QCD 



17|. 



II. OVERLAP FERMIONS COUPLED TO LARGE GAUGE THEORY 

Overlap fermions preserve chiral properties of fermions on the lattice and enables a 
study of fermionic observables in the massless limit of quarks. This makes it possible to 
compute the low energy parameters in the chiral Lagrangian of large A^ QCD. In addition, 
large A^ gauge fields on the lattice naturally come in disconnected subspaces where each 
subspace is assigned a topological charge defined using overlap fermions. 

The main point is the presence of a lattice phase transition in the lar ge N limit of pure 
gauge theory which is an extension of the Gross- Witten phase transition [20| in d = 2. The 
single plaquette Wilson action is given by 

S = bN Tr?7^,(x) (1) 

U,.{x) = U,{x)U,{x + fi)Ul{x + i))Ul{x) (2) 

where h = is kept fixed on the lattice as A^ — > oo. The continuum limit corresponds to 
b ^ oo. The lattice is taken to be a periodic torus with L sites in all directions. The link 
variables, U^{x), and the plaquette variables, U^v{x)^ are in SU(N). Let e<-(^);A;= l,---,Ar 
with — TT < 9f^i^{x) < vr denote the A^ gauge invariant eigenvalues of Ui^,y{x). Consider the 
observable, 

= I / [^^.(^)] ^^,^[^_^^ 5 " (^U^)]^'-^ Z = j [dU,{x)]e^ (3) 



5 



20 



. 6f = 0.4 
15 . It is a 



interpreted as the distribution of the plaquette eigenvalues. This distribution should be 
uniform at 6 = and should peak around 6 = in the weak coupling limit. In the weak 
coupling limit at any finite N, one expects a sharply peaked distribution at ^ = with 
p(±7r) being highly suppressed. This becomes a phase transition in the ^ oo limit at 
b = . The distribution is non-zero over the whole range of 6 if b < b^ and the distribution 
has a gap with p{6) = for \6\ > 9^{b) if 6 > 6f . The critical coupling is exactly at bc{B) for 
large L and this is a lattice transition with no continuum counterpart. The phase transition 
is at 6f = 0.5 in = 2 for all L > 1 and is the Gross- Witten phase transition 
in = 3 for all L > 3 in = 3 Q and 6f = 0.36 in = 4 for all L > 9 in = 4 
very strong first order phase transition in c? = 4 making it possible to stay in the metastable 
weak coupling side of the phase for 5 < L < 9 [iS^ . 

Gauge fields on the weak coupling side of the above lattice phase transition naturally fall 
into disconnected pieces. Since there is an upper bound on the plaquette angle, 6, there exist 
two different gauge fields, U and V say, such that it is not possible to continuously deform 
U into V without violating the upper bound somewhere in-between. One can conjecture 
that the disconnected pieces of the gauge field space as dictated by the upper bound on the 
plaquette angle are in one-to-one correspondence with the topological charge of the gauge 
field defined using overlap fermions. 

Consider a typical update of a large gauge field on the lattice using a combination 
of Cabibo-Marinari heat-bath and SU(N) over-relaxation [l^ in the weak coupling phase. 
A cold start where all gauge fields are set to unity will result in thermalized gauge field 
configurations with zero topological charge. There will be no tunneling events between 
different topological charges due to the gap in the eigenvalue distribution of the single 
plaquette. One can obtain thermalized configurations with different topological charges by 
starting from a gauge field configuration with the appropriate topological charge which also 
has a gap in the single plaquette eigenvalue distribution. For example, one could start 
with a "uniform" instanton configuration in ci = 4. The gauge field is abelian in nature 
with a uniform field strength equal to |^ in a certain plane and a certain direction in color 
space Q]. Fig. ^ shows thermalized configurations in the (5 = and Q = 1 sectors in the 
metastable weak coupling phase on a 7^ lattice with A^ = 23 at 6 = 0.35. The two histories 
correspond to two different starts where one is a cold start with (5 = and another is a 
uniform instanton start with Q = 1. No tunneling was observed in either case. The left 
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panel shows the history of the average plaquette and the right panel shows the distribution 
of the plaquette eigenvalues. There is a gap in the distribution and both of them look 
identical. It is useful to note that a naive discretization of the classical SU(2) instanton 
will not work since its plaquette eigenvalue distribution will not satisfy the gap condition. 
In practice, the energy distributions in the different topological sectors are very close to 
each other. Therefore, it is possible to offer one thermalized confieuration in one topoloeical 



0.535 




Configuration number Bin 



FIG. 1: Plaquette history (left panel) and plaquette eigenvalue distribution (right panel) in the 
weak coupling phase. Two histories, one thermalized with Q = and another with Q = 1 are 
shown at L = 7, AT = 23 and 6 = 0.35. 

The gap in the eigenvalue distribution of the plaquette has implications for the overlap 
Dirac operator in the weak coupling phase of large N QCD. The massless overlap Dirac 
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operator jl9| is defined as 



^ l + 75sign(g^(M)) 



Hw{M) is the Wilson Dirac operator at mass M, which should be chosen with —2 < M < 
and is an irrelevant parameter in the continuum limit of the theory. Naively, one expects a 
gap in the spectrum of if^(M) around zero unless M is close to the critical Wilson coupling 
and by definition one needs to use a value of M below this critical coupling above. The 
flow of eigenvalues of if^(M) as a function of M for a fixed gauge field configuration is 
directly related to gauge field topology: The net number of levels that cross zero equals 
the topological charge ^]. One can prove that a gap in the spectrum of the eigenvalue 
distribution of the plaquette results in the gap in the spectrum of if^(M) as long as one is 
significantly away from the critical coupling . Unfortunately, the eigenvalue distribution 
of the plaquette does not have a gap for = 2 or = 3 when gauge fields are generated 
using the single plaquette Wilson action. As such, the spectrum of if^(M) also does not 



have a gap below the critical coupling 



23|. But this is not the case in large QCD 



and this gives an unambiguous definition of the topological charge. Furthermore, it also 
simplifies the numerical evaluation of sign(if^(M)) which is the most time consuming part 
in a computation involving the overlap Dirac operator. 



III. PHASES OF THE LARGE N GAUGE THEORY 

As discussed in the previous section, one has to be in the weak coupling phase of large 
QCD on the lattice, defined by the presence of a gap in the eigenvalue distribution of 
the plaquette, in order to properly reproduce continuum large A^ QCD. The location of this 
lattice transition is at 6 = 6f . In addition, one has to be in a phase where none of the Zjsi 
symmetries associated with the Polyakov loops in the d directions are broken. We will refer 
to the strong coupling phase as Oh since it is a "hot" phase with all d symmetries intact. 
On the weak coupling side, there will in general he d+1 phases referred to as Oc, Ic, ■ ■ ■, dc 
where 0, 1, ■ ■ ■ , d of the d Zj^j symmetries are broken. All the phase exist as long as L > 3 in 
d=3 and L > 9 in c? = 4. Furthermore, one can be in a metastable Oc phase for 5 < L < 9 
since the lattice transition separating the Oh and Oc phase at 6 = 6f is strongly first order. 

Eguchi-Kawai reduction holds in the Oc phase and large A^ QCD is properly reproduced 
in this phase without any finite volume effects. In particular, consider a K x K Wilson 
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loop on an L'^ periodic lattice with K < L ot K > L. When K > L, the Wilson loop is 
folded since some links appear more than once in the operator. Reduction implies that the 
eigenvalue distribution of all Wilson loops are independent of L. This is illustrated jl^ in 
Fig. 121 using a 4 x 4 Wilson loop on 4'^ and 6^ lattice at a fixed coupling and fixed A^. A 
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FIG. 2: Eigenvalue density distribution of a 4 x 4 Wilson loop on 4^ (folded) and 6^ (unfolded) at 
b = 0.66 and N = 23. 



proper realization of continuum large N QCD also implies proper scaling of observables in 
the Oc phase and this is illustrated 1J| in Fig. El using a. L x L Wilson loop on an lattice 
for fixed L/b and two different L values. ' 

Unlike the Oh to Oc transition, the rest of the transitions are physical and therefore, 



the transition point will scale properly. Our numerical analysis IJ, ll5| indicates that all 
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FIG. 3: Eigenvalue density distribution of L x L Wilson loop on for L/b = 5 and L = 4, 6. 
is set to 23. 

transitions are distinct in the continuum limit but the main focus so far has been on the Oc 
to Ic transition. One useful observable to locate the phase transition is 



(5) 



where is the Polyakov loop (not its trace) in the fi direction located at a fixed point 
in the 3-plane perpendicular to yU. The e* j ; j = 1, ■ ■ ■ , are the eigenvalues of P^. 
The averaging is over all points in the 3-plane and over configurations. Equally spaced 
eigenvalues respect the Z{N) symmetry in the fi direction and maximize p to 0.5. When the 
eigenvalue spectrum starts getting modulated and opens a gap, p drops below 0.5. Fig. HI 
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shows that b = 0.351 on 6^ lattice is in the Oc phase whereas Fig. El shows that b = 0.3568 
on 7^ lattice is in the Ic phase. 

L=6 N=31 b=0.351 started from 1c 

(measurements taken every 50 lattice passes) 
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FIG. 4: History of the variable p(-P^) for each direction. We see the evolution from a state where 
one of the four Z{N) factors is broken to one in which all four are preserved. 

Let be be the location of the numerically obtained transition on an L'^ lattice. This can 
be inverted to define a critical size, Lc{b), that denotes the Oc to Ic transition. The coupling 



b has dimensions of length in d=3 and a numerical analysis [hy] shows that Lc{b) = 4.6(4)6. 
The numerically obtained |l5| critical size Lc{b) in d = A shown in Fig. [HI scales according 
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L=7 N=31 b=0.3568 started from Oc 



(measurements taken every 50 lattice passes) 
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FIG. 5: History of the variable p{P^) for each direction. We see the evolution from a state 
where all four Z{N) factors are preserved to one where one factor is broken. During the first fifty 
passes (before the first measurement) Polyakov loops in direction 3 have acquired some structure 
but, ultimately, direction 2 is selected for breakdown and the Polyakov loops in the other three 
directions converge to a symmetric state. 
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to two loop tadpole improved 2J| improved perturbation theory as 



with 



L^{b) = 0.260(15) 
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e{b) = -{TtU,Ax)) 



(6) 
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FIG. 6: The transition ranges compared to possible two loop renormalization group curves with 
tadpole improvement. 
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The Oc phase of the large continuum theory on an torus is clear and it describes 
large N QCD at zero temperature. The 4c phase can be interpreted as a continuum theory 
in finite volume and temperature with temperature equal to l/l. Rotational invariance is 
spontaneously broken in the Ic, 2c and 3c phases of the continuum theory. The physical 
interpretation of theses phases is not obvious. The Ic quite possibly describes large QCD 
in the deconfined phase since the finite size effects are only felt in the one broken direction 
and this could be interpreted as finite temperature. Then, one would conclude that large N 
QCD in the confined phase does not depend on temperature, a view that has been recently 

in a small finite box 

at zero temperature. The continuum theories described in the phases other that Oc and Ic 
might have implications not directly related to large QCD. A careful study of the Oc to Ic 
transition would tell us if the transition is the physical finite temperature phase transition. 
The critical box size associated with the Oc to Ic transition in units of string tension 
is /cv^ ~ 1.56. A comparison with the numerical estimates of the critical temperature [3| 
indicate that IrTr = 1. Recent numerical studies indicate a finite latent heat associated with 
this transition [22| providing further supporting evidence for this to be the physical finite 
temperature phase transition. 

IV. SPONTANEOUS CHIRAL SYMMETRY BREAKING 

Spontaneous chiral symmetry breaking is an important phenomenon in QCD. The mass 
of the pion made up of two light quarks is not equal or close to the sum of quark masses but 
proportional to the square root of the quark mass since the pion is a pseudo-Goldstone boson 
associated with spontaneous chiral symmetry breaking in QCD with massless quarks. Several 
lattice studies have investigated this phenomenon in the past using Wilson fermions |3| 
and staggered fermions Q] and more recently using overlap fermions . Quenched lattice 
studies of this phenomenon suffer from quenched pathologies. The quenched pathologies 
arise out of the fact that the singlet meson, r]', becomes a pseudo-Goldstone boson in the 



quenched approximation 



3l|. This introduces a new mass parameter mo in the chiral 



Lagrangian with mo ~ 900MeV in the chiral limit. The pion mass no longer is linear 
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in the square root of the hght quark mass and instead is given by 

9 ^ r ""^n 0.6 

where m is the hght quark mass and /^r = 93MeV is the pion decay constant. Spontaneous 
chiral symmetry breaking occurs due to the presence of a non-zero density of eigenvalues of 
the massless Dirac operator near zero. Let p(A) denote the density of the paired non-zero 
±zA of the anti-hermitian massless Dirac operator after removing the exact zero eigenvalues 
due to gauge field topology. The chiral condensate 83^ . 

lim lim (ipiplm)) = 7rp(0), (9) 

is non-zero if p(0) is non-zero and this shows spontaneous chiral symmetry breaking. The 



effect of the rriQ term in quenched QCD results in a diverging chiral condensate 
by 



given 



lim {'ijjip{m)) oc m (10) 

Typical quenched calculations on the lattice focus on the behavior of pion mass as a function 
of quark mass to extract the quenched divergence. Since 6 > 0, the effect of quenched 
divergence is to make the pion mass at a fixed quark mass a little heavier than what it 
should be in unquenched QCD. Since all calculations are performed at finite volume, the 
effect due to quenched divergence cannot be distinguished from finite volume effects [s^ . 
On the other hand, the chiral condensate diverges in the quenched approximation and this 
effect goes in the opposite direction from that of finite volume effects. In the large N limit, 
6 goes to zero and there are no quenched pathologies. 

In order to avoid quenched pathologies, one has to take the large limit before one 
takes the infinite volume limit. Furthermore, Eguchi-Kawai reduction holds in the Oc phase 
and therefore the infinite N limit in this phase will not depend on the physical volume. We 
should therefore see evidence for a chiral condensate on a finite lattice in the large N limit 
as long as we are in the Oc phase. This is counter-intuitive since we expect to see a finite 
density of eigenvalues near zero in the spectrum of the massless Dirac operator even at finite 
lattice volume. Chiral random matrix theory jsJ] provides an understanding of spontaneous 
chiral symmetry breaking in large QCD on a finite lattice |16i] . 

Simple counting of degrees of freedom shows that (ipip) is proportional to A^. Therefore 
the low lying eigenvalues of the massless Dirac operator should scale like if chiral sym- 
metry is spontaneously broken in the large A^ limit. This is a consequence of level repulsion. 
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The paired non-zero eigenvalues, ±A, of the massless hermitian overlap Dirac operator sat- 
isfies < A < 1. The operator has 2NV positive eigenvalues (assuming zero topology) on a 

V = lattice. Level repulsion would result in eigenvalues that are roughly equally spaced. 
Therefore, the mean spacing will be proportional to and the lowest eigenvalues will scale 
like j^. Furthermore, the scaled eigenvalues z = ASA^V^ should essentially be distributed 
in some universal manner and this universal behavior is dictated by chiral random matrix 
theo.. 0. The .ale, .ha. connec. . a„d A . the cM.a, conde.a.e. Unive.a, beha.o. 
is approached as NV, the number of eigenvalues, go to infinity. This can be achieved by 
keeping fixed and taking V to infinity which is the conventional approach. But, it can 
also be achieved by taking to infinity at fixed V and this makes use of the reduction 
in large QCD. Quenched divergences imply that A does not scale like at fixed A^ as 

V oo. In spite of this, it has been shown that one can define a scale ^{V) such that chiral 
random matrix theory is respected in the quenched approximation 3^. But, S(V^) diverges 
as V ^ oo. 

One could show evidence for spontaneous chiral symmetry breaking in large A^ QCD 
by a direct computation of S in the Oc phase. This can be done by a computation of 
lim^^o liKiiv^oo(V'V'(^)) oil ^ finite L'^ lattice. In this case, it is necessary to take the large 
A^ limit before taking the massless limit since the chiral condensate is zero at m = for 
any finite A^. Numerically, one has to perform some sort of a finite A^ analysis equivalent 
to the usual finite volume analysis done in order to extract an order parameter in the 
thermodynamic limit. Chiral random matrix theory provides us with a way to perform the 
finite A^ analysis without a computation of the chiral condensate. Let Aj, i = 1, ■ ■ ■ , 2NV, 
with Ai < A2 < ■ ■ • A2Ary be all the positive eigenvalues of the hermitian overlap Dirac 
operator. Let = Xk^NV be the scaled eigenvalues that obey the universal distribution 
given by chiral random matrix theory js^. As A^ increases, more and more zj. will obey 
the universal distributions. In particular, one will find that p{r = ^) obeys a universal 
distribution as long as A^ is large enough. Once this is achieved, S can be obtained by 
matching (zi) and (22) to their values given by chiral random matrix theory. A numerical 
computation shows that one can obtain S in the large A^ limit of QCD using chiral random 
matrix theory jisj. Existence of a chiral condensate in the continuum limit amounts to 
showing that obtained on the lattice scales properly in the continuum limit and one finds 
that EV3/^ = 0.588 which translates to S^/^ ~ l55MeV if we assume that f = = 264MeV. 
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A. Chiral condensate in the large N limit of QCD2 



Before proceeding to the computation of chiral condensate using chiral random matrix 
theory, it is useful to verify the ideas of chiral random matrix theory in two dimensional 
QCD jl6|. The chiral condensate in the 't Hooft model is exactly known j38|| and it is given 
by 

m - ^ (u) 

VOTTO 

for the overlap Dirac operator, where M is the Wilson mass parameter used in the definition 
of the overlap Dirac operator. Eguchi-Kawai reduction holds for all L in c? = 2 and therefore 
one can compute the chiral condensate by taking the large limit on a 1^ lattice. One 
can compute the full spectrum of the overlap Dirac operator for very large N using exact 
diagonalization on a 1^ lattice. Let pi{zi) and P2{z2) be the universal distributions of the 
two lowest eigenvalues and let p{r) be the distribution of the ratio of the first eigenvalue and 
the second eigenvalue. Fig. [7| shows that approach to chiral random matrix theory as is 
increased on a 1^ lattice at a fixed b. N = 37 is too small for chiral random matrix theory to 
be satisfied. At = 79 the first eigenvalue already satisfies the universal distribution and 
both eigenvalues satisfy the universal distribution beyond A^ = 137. The coupling was set 
to 6 = 1 and we used M = — 1 in the definition of the overlap Dirac operator. Therefore, we 
should find S = 0.23 as A^ ^ 00. Furthermore, < r >= 0.37 as A^ — > 00. The approach to 
the N —>■ 00 limit of < r > and S using the first and second eigenvalue is shown in Fig. |H1 
One can see that the first eigenvalue approaches chiral random matrix theory before the 
second eigenvalue. 



B. Chiral condensate in the large N limit of QCD 

Analysis of QCD2 shows that one can use chiral random matrix theory to extract the 
chiral condensate in the large A^ limit of QCD. To do this, one has to be in the Oc phase of 
QCD. Therefore, we worked with L in the range of 6 to 10. At each L, there is a range of 
coupling < b < bc{L) where the theory is in the Oc phase. It is best to pick the coupling 
close to bc{L) to minimize lattice spacing effects at that L. We therefore used b = 0.350 
for L > 6, 6 = 0.355 for for L > 7 and b = 0.3585 for L > 8. We worked with L = 6, 7, 8 at 
b = 0.350 to ensure that reduction holds. We worked with L = 8, 9, 10 at 6 = 0.355 again 
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FIG. 7: Approach to chiral random matrix theory on a 1^ lattice at 6 = 1 as is increased from 
37 to 273. 
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FIG. 8: Estimation of the chiral condensate in the N ^ oo hmit of QCD2 at 6 = 1 as A'^ is increased 
from 23 to 273. 

to ensure that reduction holds. We worked with L = 9 at 6 = 0.3585 and we also obtained 
one result on a coarse lattice of 6 = 0.346 at L = 9. The scaling properties of the code are 
as follows. The gauge field generation scales like N^L^. The computation of the lowest two 
eigenvalues of the hermitian overlap Dirac operator scales like N^L^ and this arises from 
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two factors. One factor is due to the action of the Wilson Dirac operator on a vector and 
this scales as N'^L^. The low lying eigenvalues of the overlap Dirac operator are obtained 
using Ritz techniques and the relevant condition number is the ratio of the highest to the 
lowest eigenvalue and therefore the second factor scales as NL'^. 

Fig. El shows the approach to chiral random matrix theory as one approaches the large 
N limit at a fixed L and b. One sees a clear agreement with chiral random matrix theory in 
the distribution of the ratio of the two lowest eigenvalues for > 23. All configurations are 
in the zero topological sector. The universal distribution of p{r) is different in the different 
topological sectors and Fig. El shows that one finds agreement with chiral random matrix 
theory in the Q = 1 topological sector when there is agreement in the Q = topological 
sector. Finally, Fig. ^2 shows that one also finds agreement with chiral random matrix 
theory as one goes toward the continuum limit. 

The first and second eigenvalues are used to obtain the chiral condensate, E, by setting 
the average of Zi = Xi'EiNL'^; 2 = 1,2 equal to the ones dictated by chiral random matrix 
theory. One should find Si = S2 for values of b, L and where p{r) agreed with the 
universal distribution given by chiral random matrix theory. Fig. El shows the results for 
El and S2 from all our numerical simulations at 6 = 0.350 where we found p{r) to agree 
with chiral random matrix theory. First of all we find that Si = S2 for all the data shown 
in Fig. El Furthermore, we do not see any finite volume effects at 6 = 0.350 indicating that 
reduction works properly for fermionic observables. The finite A^ effects are relatively small 
over the range of A^ shown since the y-axis in Fig. IT^ shows a deviation of only 5%. We can 
therefore estimate S^/3(6 = 0.350) = 0.142 ± 0.006. 

Estimation of the chiral condensate using chiral random matrix theory has proved to be 

n 

far superior to conventional techniques [IQ] . The main problem that plagues the conventional 
approach is the quadratic divergence in (ipip). Therefore, large subtractions are needed in 
order to extract the chiral condensate. This can remedied by considering derivatives with 
respect to of the chiral condensate. This ameliorates the problem but the dependence 



on small quark masses was stil 
condensate using extrapolation 



strong and this makes it difficult to estimate the chiral 



16|. In addition to the need for subtractions, one can only 



obtain a stochastic estimate of this quantity and it is necessary for the random source to 
have good overlap with the low lying eigenvalues. On the other hand, numerical evaluation 
using Ritz functional techniques results in estimates of low lying eigenvalues with errors that 
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FIG. 9: Approach to chiral random matrix theory at a fixed b and L with increasing N. 



are less than a tenth of a percent. 

If we interpret the Oc to Ic transition as a finite temperature phase transition, then one 
would expect chiral symmetry to be restored in the Ic phase. But, chiral random matrix 
theory arguments would indicate that chiral symmetry would still be broken since the 
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FIG. 10: Evidence for agreement with chiral random matrix theory in the Q = and Q = 1 
topological sectors at a fixed b = 0.350. 



22 



3.5 



3 — 
2.5 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

L=8, N=23, b=0.350, 257 cfgs — 



3.5 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

3 I— L=8, N=23, b=0.355, 288 cfgs - 




_ L=9, N=17, b=0.355, 288 cfgs _ 



3.5 
3 

2.5 I — 



0.2 0.4 0.6 0.8 

r 

I I I I I I I I I I I I I I I I I I I I I I I I I I I I 




3.5 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

L=9, N=17, b=0.3585, 336 cfgs — 




3.5 

3 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

- L=10, N=13, b=0.355, 288 cfgs 




FIG. 11: Evidence for agreement with chiral random matrix theory as one approaches the contin- 
uum limit. 
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FIG. 12: Estimate of the chiral condensate in the large N Umit of QCD. 



symmetry is broken only in one of the four directions. Three other directions still behave as 
if they were infinite in extent and one would still expect level repulsion with level spacings 
that go like Therefore, one would expect chiral symmetry to be broken in the Ic phase 
with a value that depends on the length of the direction which is broken. Therefore, it 
would be interesting to compute the chiral condensate in the Ic phase. Similar arguments 
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would also say that chiral symmetry is broken in the 2c and 3c phase. Chiral symmetry 
is expected to be restored only in the 4c phase. A study of chiral symmetry breaking in 
these other phases would be interesting on their own right for a full understanding of the 
associated continuum theories. 



V. CONCLUSIONS AND FUTURE WORK 



Phase transitions as a function of scale are ubiquitous in large QCD. Continuum 
reduction says that large A^ gauge theory on an l'^ torus undergoes a phase transition 
when I = Ic and large A^ QCD is properly reproduced as long as / > Ic- This can used to 
numerically solve large A^ QCD by working on small lattices. Furthermore, one can work in 
the quenched approximation on the lattice. It is important to take the N —>■ oo limit before 
one takes any other limit. Then, one finds that there are no finite volumes effects and it is 
possible to obtain the low energy parameters of the chiral Lagrangian in the large A^ limit of 
QCD 3]. Chiral random matrix theory proved to be a useful tool in the extraction of 
the chiral condensate since it describes the behavior of small eigenvalues for large but finite 
number of degrees of freedom. 

The next step is the computation of the pion decay constant and work is close to com- 
pletion This will establish the validity of momentum quenching in the computation of 
meson propagators. A study of current correlators will be useful in the extraction of the 
mass of the p. 

On the technical side, the Wilson Dirac operator has a gap in its spectrum at values of 
M used in the overlap Dirac operator. This is due to the lattice phase transition present in 
the large A^ limit of QCD and is a consequence of the gap in the eigenvalue distribution of 
the single plaquette. Therefore, gauge fields come in disconnected subspaces. In addition, 
the gap in H^iM) results in a significant reduction in the computational cost of the action 
of the overlap Dirac operator on a vector. 

The Ic to 4c phases discussed in section 3 have possible implications for certain string 
theories ^ and it would be useful to compute the chiral condensate in all these phases. 
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